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We study the thermodynamics of QCD in the limit of large flavor number (Nf) and test 
the proposal to resum the physics of hard thermal loops (HTL) through a nonperturba- 
Pnj tivc expression for the entropy obtained from a ^-derivable two-loop approximation. The 

^ ■ fcrmionic contribution to the entropy involves a full next-to-leading order evaluation of the 

asymptotic thermal quark mass, which is non-local, and for which only a weighted average 
^^"O ■ value was known previously. For a natural choice of renormalization scale we find remarkably 

good agreement of the next-to-leading-order HTL results for the fcrmion self energy and in 
turn for the entropy with the respective exact large- Nf results even at very large coupling. 

(N : 

m ■ 

O ' I. INTRODUCTION 

o\ : 
o . 

!/"■) . The perturbative series for the thermodynamic potentials of hot QCD is by now known up to 

and including order g 6 log(g) 0, 0, 0, 0] , with g = ^4ira s being the Yang-Mills coupling constant. 
Taken at face value, this series is poorly convergent and suffers from a strong dependence on 
the renormalization point even at temperatures many orders of magnitude higher than the QCD 
Ph! scale Aqcd- However, this problem is not specific to QCD at high temperature, which has a 

nonperturbative sector starting to contribute at order g®. Similarly poor convergence behaviour 
appears also in simple scalar field theory 0], and even in the case of large- N ^ 4 theory jfj, where 
all interactions can be resummed in a local thermal mass term, as soon as one starts expanding 
^ ' out in a series of powers and logarithms of the coupling. 

■ Various mathematical extrapolation techniques have been tried to restore the convergence of 

the perturbative series, such as Pade approximants 0,0,0, self-similar approximants [u|, and 
Borel resummation A more physically motivated proposal for reorganizing thermal pertur- 

bation theories is called "screened perturbation theory" [13. This is a variant of a variational 
perturbation theory, where the tree-level Lagrangian is modified so that it includes a mass term, 
which is then determined by a variational principle (minimal sensitivity) . In the case of gauge the- 
orieSjthe prefactor of the gauge-invariant nonlocal and nonlinear hard-thermal-loop (HTL) action 
llq ] is used for this purpose in a generalization of this approach to QCD by Andersen, Braaten, 



Ph! 
P. 



Petitgirard, and Strickland [12 



m a generali 



The HTL action is the correct leading-order effective action for soft modes at energies para- 
metrically smaller than the temperature. However it is somewhat problematic to also use it at 
hard scales, which are responsible for the leading terms in the thermodynamic potential, albeit the 
problems with thermal perturbation theory clearly come from screening effects at soft scales. In- 
deed (HTL-)screened perturbation theory changes the ultraviolet structure and requires an ad-hoc 
renormalization of additional UV singularities. 

An alternative proposal for resumming the physics of hard thermal loops was put forward in 
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Refs. 0, [H Q 1 and is based on a nonperturbative expression for the entropy density that 



can be obtained from a ^-derivable two-loop approximation |2fj]. Here the emphasis is fully on 
a quasiparticle picture, whose residual interactions are assumed to be weak after the bulk of the 
interaction effects have been incorporated in the spectral data of the quasiparticles. As opposed 
to methods based on "screened perturbation theory", ^-derivable approximations have better 
properties with respect to renormalization. They are renormalizable in the case of scalar field 
theories |27l Eil . E^ . 13(1 Iflll > whereas in gauge theories the question of renormalizability remains 
still opened. However, when combined with HTL approximations, they define UV finite physical 
quantities without introducing spurious counterterms 21, 2^, 2^ |. 



These approaches indeed succeed in taming the plasmon term ~ g 3 that spoils the apparent 
convergence of strict perturbative expansions in g and comparison with existing lattice data suggest 
that the entropy of QCD for T > 3T C can indeed be accounted for remarkably well when the 
leading-order interactions are resummed into spectral properties of HTL quasiparticles 2jl 2^, 2^ |. 



Since this success is gauged mostly from the comparison with lattice data, it is desirable to 
have other possibilities for testing these resummation prescriptions. In Ref. [331 ] it was suggested 
to use the large flavor- number (Nf) limit of QED and QCD for that purpose. The thermodynamic 
potential has been worked out in Ref. 0, to order JVj and to all orders in g Nf, which is 
kept finite as the limit Nf — > oo, g — ► is taken. 

Large- Nf QCD is no longer asymptotically free. In this limit the renormalization scale depen- 
dence is determined (non-perturbatively) by the one-loop beta function according to 



9eM 6?r 

where we have defined g 2 s = g 2 Nf/2 (in QED we would have g 2 s = g 2 Nf). There is a Landau 
singularity in the vacuum gauge field propagator at Q 2 = with 

A L = £Mse 5/6 e 67r2/ <& (/iMs) , (2) 

where /%jg is the scale of modified minimal subtraction. However, although the large- Nf theory has 
an unavoidable ambiguity associated with its UV completion, in the thermodynamic potential the 
magnitude of the ambiguity is suppressed by a factor (T/Al) 4 , which in practice means that the 
Landau pole problem can be ignored for all coupling g c s(jj.y^ = vrT) < 6. Since strict perturbation 
theory appears to work only for g e g < 2, this gives enough room for testing improvements of 
thermal perturbation theory. 

In this paper we shall test the HTL resummation proposed in Refs. 2^, [3] which use the 



nonperturbative expression for the entropy obtained from a ^-derivable two-loop approximation. 
In the large- Nf limit, a complete HTL resummation involves in particular the evaluation of next- 
to-leading-order thermal quark masses at asymptotic hard scales. In one-loop HTL-resummed 
perturbation theory, these thermal masses are nonlocal, i.e. functions of momentum, and so far only 
a certain weighted average was known, and this only to next-to-leading order in the coupling. In this 
paper we shall perform a complete HTL evaluation, which resums an infinite series in the coupling, 
and we compare with the numerical evaluation of the exact large-Ay values. A corresponding 
evaluation of the entropy shows considerable improvement compared to a previous evaluation using 
the averaged asymptotic masses at next-to-leading order in HTL-resummed perturbation theory 
0| . Although at larger coupling the HTL results have a large renormalization scale dependence, 
choosing the scale of fastest apparent convergence re pro duces the exact large- Nf result for all g 2 s 
with a quality comparable to an optimised g 6 result |3J|. 



1 For an attempt to achieve a comparable resummation directly in terms of the pressure see $Z 
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The organization of this paper is as follows: In section[H]we first review the ^-derivable approx- 
imation to the entropy. At two-loop order we present the remarkably simple quasiparticle formulae 
for the entropy and number density which can be used both in the large- Nf limit and for an HTL 
approximation. In section II I II we discuss the issue of renormalization both from the perspective 
of the thermodynamic potential and of the quasiparticle entropy and number density. In section 
IIVI we evaluate the fermion self-energy, which plays a central role in the quasiparticle entropy ex- 
pressions, and which requires to solve a technical problem related to the necessity of eliminating 
the Landau pole by a momentum cutoff. Such a cutoff has to be imposed in a way which respects 
Euclidean rotational invariance in order to avoid spurious contributions in already renormalized 
expressions [derived in dimensional regularization as we shall use modified minimal subtraction 
(MS)]. Since the quasiparticle entropy is formulated in Minkowski space, this requirement leads 
to intricacies in the numerical evaluation, whose results are finally given in section There we 
present the numerical results for the complete large-iVy entropy and the complete momentum- 
dependent asymptotic thermal quark masses together with their respective HTL approximations. 
We find remarkable agreement of the latter with the exact results provided the HTL approximation 
is used in a nonperturbative manner (i.e., not truncating at the order of perturbative accuracy) 
and an optimized renormalization scale is chosen. In this case, the agreement is comparable to 
optimized results using perturbative results through order g e . 

At this point it is useful to specify the notation to be used throughout the paper. As gauge 
group we take SU(iV), so N denotes the number of quark colors. As already said, Nf is the number 
of quark flavors. The number of gluons is denoted by N g = N 2 — 1. We use the Minkowski metric 
guv _ ^_|_^ _^ _^ anc [ denote by capital letters 4-momenta P with components P^ = (po,p), 
so that P 2 = p§ — p 2 . In finite temperature calculations, we are led to set po = iu> n , where 
oj n is a Matsubara frequency (u n = nir[3 where n is an even integer for bosons and odd integer 
for fermions). We shall denote the typical sum- integrals that occur at finite temperature by the 
following shorthand notation: 

At zero temperature, the sum over Matsubara frequencies is replaced by an integral over the 
imaginary energy axis, leading to an Euclidean integral denoted by 

For an integral over the four components of a four-vector, we simply write 



dpo f d 3 p 



We denote the free propagator by Aq(P) = —1/(P 2 — M 2 ) or, in its spectral representation, by 

A (u,p)= [ (6) 
Jpo Po ~ u 

where lo is a complex energy variable, and the spectral function is po(P) = 27re(j>o)<5(-f >2 — M 2 ). 
The free fermion propagator is then given by Sq{P) = {Jf — M)Aq(P) or, in its spectral form, by 

S {uj,p)= / , (7) 

Jpo Po ~ w 

where we used the fact that po(P) is an odd function in po to replace ujjq by po7o i n within the 
integral. 
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FIG. 1: Diagrams for <I> at 2-loop order in QCD. Wiggly, plain, and dotted lines refer respectively to gluons, 
quarks, and ghosts. Only diagram (d) contributes at next-todeading order in the large-iVj limit. 



II. 4>-DERIVABLE TWO-LOOP ENTROPY AND HTL RESUMMATION 



A. Generalities 



^-derivable app roximations [3?1 ] are constructed from the two-particle-irreducible (2PI) skeleton 
expansion jjjl . l39§ | . In the latter, the thermodynamic potential is expressed in terms of dressed 
propagators (G for bosons, S for fermions) according to 

n[G,S] = iTTr logG^ 1 - iTTrnG 

-TTr logS" 1 + TTrS5 + r$[G,5] , (8) 

where "Tr" refers to full functional traces Tr = tr /(f dr J x -» (3V tr f K and S] is the sum of 
2-particle-irreducible "skeleton" diagrams. In gauge theories one either has to assume a ghost-free 
gauge such as temporal axial gauge or to include the bosonic Faddeev-Popov ghost propagator 
G g h- For simplicity we shall assume a ghost-free gauge in the general discussion. At both the level 
of HTL approximations and in the large- Nf li mit, we can and shall use the Coulomb gauge for this 
purpose instead of the somewhat problematic [4(1 ] temporal axial gauge. 
Standard contour integration gives 

tt[G,S]/V = tr / n(k ) Im \ \og G" 1 (K) -U(K)G{K) 
Jk l j 

+ 2tr / f (k ) lm\ log S^i^-^^SiK)] +T$[G,S\/V, (9) 

where "tr" denotes a trace over discrete labels only. The Bose-Einstein and Fermi-Dirac factors 
are defined as n(ui) = {e u) l T — and /(a;) = (e^~^^ T + l) -1 , respectively. 

The self-energies II = G" 1 — Gq 1 and S = S~ 1 — Sq 1 , where Go and Sq are bare propagators, are 
themselves functionals of the full propagators. They are determined by the stationarity property: 

sn[G,s]/SG = o = sn[G,s\/ss , (10) 

according to 

6$[G,S]/6G = |n, 5$[G,S]/6S = — E. (11) 

The ^-derivable two-loop approximation consists of keeping only the two-loop skeleton diagrams 
(see Fig. , which leads to a dressed one- loop approximation for the self-energies (fTTf) . In a gauge 
theory this generally introduces gauge dependences (which are parametrically suppressed, though 
|4lj). However the further approximations put forward in Refs. [2ll. E3. E3| are manifestly gauge 
independent in the ^-derivable two-loop approximation. 
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A self-consistent two-loop approximation for $7 has a remarkable consequence for the first deriva- 
tives of the thermodynamic potential, the entropy and the number densities: 



c_ dP 
6 " dT 



dP 



P 



-n/v. 



(12) 



T 



Because of the stationarity property (|1()[) . one can ignore the T and /x dependences implicit in the 
spectral densities of the full propagators, and differentiate exclusively the statistical distribution 
functions n and / in Q. Now the derivative of the two-loop functional T&[G, S] at fixed spectral 
densities of the propagators G and S turns out to cancel part of the terms Im (JIG) and Im (SS) 



m 



3(T$ 2 



2— loop 



-loop J 



dT 



G,S 



+ tr 



A' 



dT dT 



Re £ Im S 



(13) 



leading to the remarkably simple formulae |22j, 12.1 12_ 



A' 



-tr 

-2tr 

-2tr 



dn(fco) 
K dT 
df(ko 



k dT 
df(ko) 



K 



Im logG" 1 - ImllReG 
Im \ogS~ 1 - ImSReS" 
Im log5 _1 - ImSReS 



(14) 
(15) 



where we have dropped the label "2-loop" that could be attached to S and M. 

Through these formulae, effectively one- loop integrals, all interactions below perturbative order 
g 4 are summarized by spectral data only, which shows that entropy and density are the preferred 
quantities for a quasiparticle description. In particular the term g 3 , which usually spoils the 
apparent convergence of strict perturbative expansions, is now incorporated in a nonpolynomial 
expression, together with (incompletely resummed) higher-order terms. 

Moreover, these expressions are UV finite as soon as the self-energies are, and thus the former 
are useful as a starting point for further approximations. In Refs. 0, 0, E^l it was proposed to 
use the gauge-invariant hard thermal loops for this purpose. Before considering HTL's, we recall 
some general properties of self-energies. 



B. Bosonic and fermionic self-energies 

In the Coulomb gauge, the gauge propagator can be decomposed into a longitudinal and a 
transverse contribution 

Gqo(^) = G L (K), 

G^(K) = { Sij -^}G T (K), (16) 
and we define corresponding self-energy components through 

a ^ K) = wrhK)- CT < K >- -Wn T( Ar (17) 



Longitudinal and transverse spectral functions are introduced according to 

G^,k) = -^ 2+ f f®, G T (.,k)=f^l. (18) 
k 2 Jk a k -uj J ko k -uj 

In the following it will be convenient to write the temporal component of the Coulomb gauge 
propagator alternatively as 

Goo(K) = H^.G e (K), (19) 

and the Dyson equations as 

g^ 1 = -k\ + k 2 + n t ,th + n iiVac , 

Gt = ~k 2 o + k 2 + U £)th + Ue tVac , (20) 

such that IIt = n t th + n t)Vac but n^th + n^,vac = — ^(fcg — k 2 )/k 2 . Here we have separated off 
the vacuum (T, /z — » 0) limit of the self-energy components. In the large- Nf limit, we shall have 
that IIt )Vac = Il^vac = n vac even in the non-covariant Coulomb gauge. 

The most general form of the self-energy £ at finite temperature and density can be written as 

E(K) = a(iY) 7 ° + b(K)k • 7 + c(K) . (21) 

We define the projection of the self-energy on — M as 

t{K) = tr [(#' - M)Y,{K)} = 4wa(u, fc) + Akb(uo, k) - AMc(lo, k) . (22) 

In the massless case, c{K) = 0, and the quark propagator at finite temperature or density can be 
split into two separate components with opposite ratio of chirality over helicity: 

^ ± i*U(ior (23) 

where 

E±(K) = b(K) ± a(K) . (24) 
Furthermore, S(if) and E-^if) on the light cone uj = ±k are related by 

S±(w = ±fc, k) = ^E(w = ±ife, k) (25) 
(see also Appendix El for a relation between S + and £_). 

C. HTL approximation of the self-consistent entropy 

The HTL effective action [lil Il6| is an effective action for soft modes with energy scales ~ 
gT, which are, at least parametrically, at smaller energy than the hard modes defined by the 
temperature scale T (or the chemical potential when /i> T). At hard scales, the HTL effective 
action is no longer accurate, except at small virtuality, but this is indeed the phase space domain 
which contributes the leading-order interaction terms oc g 2 in the expressions (|14|) and (|15j) . The 
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order g 2 contribution is obtained by expanding out the propagators and keeping a single self-energy 
insertion. In SU(iV) gauge theory with N g = N 2 — 1 gluons and Nf massless quarks this leads to 



2N„ 



dn(k ) 



+4NNf 



dT 

df(k 



Re IIt Im 



h 2 — h 2 

rt'Q r\j 



dT 



Re Im 



kn — k 



Re £_ Im 



1 



k + k 



(26) 



and similarly in the density expression. The imaginary part of the free propagator puts the self- 
energy insertions on the light-cone, where they are accurately (to order g 2 ) given by the HTL 
value 2 of the transverse component of the gluon self-energy IIt 



U T (k = k) 



\rh 2 D , 



(27) 



and the HTL value of £± = Mjk/\k\ ± 7o)S, 



£±(k = ±k) 



M 2 
2k 



~k~ 



(28) 



even though k is no longer soft |42], |43j • The leading-order terms in the interaction contribution to 
the thermodynamic potentials S and N are thus related to the "asymptotic" thermal masses jtt-oo 
and Mqo of hard gluons and fermions, respectively. As indicated in Eqs. (f2*7|) and (|2*8|) . at leading 
order these are proportional to the HTL Debye mass mo and HTL fermionic plasma frequency M, 
given by 



m 2 D 



2tt 2 



M 



9 2 C 



f 



T 2 + 



7T- 



(29) 



for (uniform) quark chemical potential /i and Cf = (N 2 — 1)/(2N). 

However, higher-order contributions in the thermodynamic potential cannot be calculated by 
expanding out the self-energies in the expressions for the thermodynamic potentials given in Section 
III Al as this would lead to infrared divergences in the electrostatic sector. At least the Debye 
mass that appears in the static propagator needs to be resummed, after which one can calculate 
systematically up to and including order g e log g. But already the first contribution from the soft 
sector, the "plasmon term" oc g 3 spoils the apparent convergence of a strict perturbative series. 

With the above nonperturbative expressions for entropy and density it is now possible to resum 
simultaneously the effects of Debye screening and other collective phenomena, such as dynamical 
screening and Landau damping. At soft momentum scales, these are determined to leading order 
by the HTL self energy expressions 



and 



n L {k ,k) 
U T (k ,k) 

£±(k ,k) = 



m 2 D 



1 



A^-lo fco + k 
2k ko — k 



k 2 

m D -\ 



k 2 



k 2 



n L 



M 2 
~k~ 



ko^k k 
1 — log 



2k 



kn — k 



(30) 



(31) 



2 Here and in what follows we shall denote HTL values by putting a hat on the quantities in question. 
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FIG. 2: NLO contributions to STIt at hard momentum. Thick dashed and wiggly lines with a blob represent 
HTL-rcsummed longitudinal and transverse propagators, respectively. In the large-iV/ limit the blobs rep- 
resent full 1-loop resummcd propagators. In this limit, these diagrams are suppressed by a factor of 1/Nf. 



FIG. 3: Same as Fig. [3] for the NLO contributions to ST, at hard momentum. Although in the large- Nf 
limit these diagrams are also suppressed individually by a factor of Nf, they nevertheless contribute to the 
entropy at order N®. 



Now Dyson-resumming these self-energies into dressed propagators and inserting them into 
the entropy formula (|14|) turns out to account for only a fraction of the plasmon term ~ g 3 . In 
the entropy formula (|14[) the larger part of the plasmon term arises from hard momentum scales, 
namely from corrections to the leading-order asymptotic masses. 

Calculation of the next-to-leading-order corrections to the asymptotic thermal masses (|27|) and 
(|28j) itself requires HTL resummation. However, because they are given by self-energies with a 
hard external momentum, they involve only a single HTL propagator and no HTL vertices, see 
Figs. Hand El 

These corrections to the asymptotic thermal masses are, in contrast to their leading-order (HTL) 
values, nontrivial functions of the momentum, 

Srn^ik) = Re6U T (k = k ), $ M L( k ) = Re2k5Z + (k = k). (32) 

At next-to-leading order HTL perturbation theory these expressions can be evaluated only numer- 
ically. However, through their contribution to the plasmon term, the following weighted averages 
are determined to order g s and given by remarkably simple results 21, 2^, 2^] 



fdkk^l " 2tT 



5m oo = ~ , fbgfcF = —fNTm D + 0( 5 4 ) , (33) 

ar' 

- „ f dkk^-8WL(k) 1 , , ^ 

smi = j I ; k T kd _nf [) = + o{9) • (34) 

Note that these results pertain only to the hard excitations; corrections to the various thermal 
masses of soft excitations are known to differ substantially from (|33|) . For instance, the relative 
correction to the gluonic plasma frequency 01 at = 0, <5m^ /m^ , is only about a third of 
<5 m oo/ m oo) the NLO correction to the nonabelian Debye mass on the other hand is even positive 
for small coupling and moreover logarithmically enhanced [iil ]. 

Sm 2 D = +^-g 2 NTm D log -, (35) 
2tt g 

where the constant under the logarithm is nonperturbative and cannot be calculated by weak- 
coupling techniques. However these corrections to the dispersion laws at soft momenta lead to 
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contributions to the thermodynamic potential which are beyond the perturbative accuracy of the 
two-loop functionals (|14j) and (|15j) . and thus we do not include them in our definition of next-to- 
leading-order HTL approximation of the entropy or density. This we instead define as an evaluation 
of p4jl and l|15|) with HTL propagators and self-energies, which at hard momenta include the 



complete corrections to the asymptotic thermal masses (|32j) . In Refs. 2l|, l22l . 12.1 l3q only the 
averaged next-to-leading order asymptotic masses (j33[l and (j34[) were taken into account. 

In this paper we shall evaluate the complete momentum dependence of the next-to-leading- 
order asymptotic quark masses <5M^(k), which is the only correction also relevant at large Nf. 
While both (|33j) and (|3"1|) are suppressed by 1/Nf in the large- Nf limit, the latter eventually gets 
multiplied by a factor Nf. We then compare with the exact nonperturbative results that one can 
derive in the large- Nf limit 



III. LARGE- iV f EXPANSION OF THE THERMODYNAMIC POTENTIAL 



We shall now focus on the large- Nf limit of the thermodynamic potential, which is nonpertur- 
bative in the effective coupling constant and therefore requires nonperturbative renormalization. In 
this respect we extend the previous derivations of Ref. to include nonzero quark masses 

(see also |!f3]), even though the final numerical evaluation will be carried out only for the massless 
case. We then discuss the special features of the self-consistent entropy and density expressions in 
the large- iV^ limit. 



A. Large- Nf limit of the 2PI thermodynamic potential 



The thermodynamic potential of Eq. (jHJ) is given in terms of full propagators G and S and 
the <3? functional, the sum of all 2PI diagrams. In the large-iVj- limit, the $ functional reduces to 
a single two-loop skeleton diagram, namely the last one of those displayed in Fig. Q Only this 
diagram remains of order 1 when g 2 — > 0, Nf — > oo with fixed g 2 s oc g 2 Nf ~ 0(1)- Furthermore, 
the fermion propagator in this diagram can be replaced by the free one, since a fermion self energy 
is of order £ ~ 1/Nf. Hence 

$[G,5] -> $*[G] = iTrOFG) (36) 

where n* is a fermion bubble involving the free fermion propagator So, and is proportional to g 2 s . 
From a variational point of view, n* is now a fixed quantity, and the self-consistent n defined 
through Eq. ((TTj) is equal to it because 3>*[G] is linear in G. Hence, n = 25Q*[G]/5G = n*. On 
the other hand, S = Sq is no longer a variational quantity at all. We thus have 

U[G,S] = -TTrlogS p ; 1 + iTTrlogG- 1 -iTTr(nG) + ^[G], (37) 

which, when evaluated at its extremal point and setting G" 1 = Gq 1 + n*, reduces to 3 

n*[G*] = -PV = — TTr IoeSq 1 + |TTr logG^ 1 . (38) 

Diagrammatically, the pressure at leading order (Nj) is thus just given by the free fermion loop. 
It contains an overall, temperature independent divergence which can be eliminated by subtracting 



3 We stress here that viewed as a functional of G* — > G, the right hand side of Eq. 138H is not stationary at 
G^ 1 = Gq 1 +11* and thus one cannot use the same simplifications to compute the entropy as in the 2PI framework. 
In order to do so, one should rather work with Eq. I|37|l 
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FIG. 4: Fcrmion bubble resummation contributing to the pressure to next-to-leading order in the 1/Nf 
expansion. 



the corresponding vacuum pressure: 



pj = NN f tr 



log 5, 



-i 



A" 



log So" 1 



(39) 



In Eq. 1)39)) and in the rest of this section 1)111 A|) . we perform explicitly flavor and color traces. Thus 
the symbol "tr" only denotes a trace over spin indices. The next-to-leading order contribution 
~ N® is obtained from the resummation of the bosonic self-energy II (see Fig. on the gluon 
propagator (ring resummation): 



nf (Q) = 9 c 2 ff tr £ 7 M S (K) ^ S (K - Q) , 
which one can split into a vacuum and a thermal piece (see Appendix I A 1)1 : 



(40) 



ac (Q) = g 2 s tr f ^ S (K) Y S (K - Q) , 
n t T(Q) = 2g 2 cS tr / ^ " M) a (K) ^ S (K - Q) , 

JK 



(41) 



With 0-o(if)=€(fco)/(|fco|)po(^). 

This function &o{K) is a particular combination of the thermal factor and the spectral de nsity 
which appears systematically when isolating thermal contributions in the Matsubara formalism (47 ] . 
The thermal factor /(|A;o|) only involves positive energies and thus the thermal contribution is UV 
finite. We will also need the dominant UV asymptotic behaviour of n^. Since the momentum K 
is cut off by the temperature, the leading asymptotic behaviour at large Q is a priori determined by 
So(K — Q) ~ l/Q- However, since the gluon self-energy is an even function of Q, this asymptotic 
behaviour is in fact improved to l/Q 2 . 4 

The vacuum piece n^ vac is UV divergent (hence the label "b" standing for "bare"). However, 
to the same order in Nf, one can add the counterterm 5Za{Q 2 — Q^Q"}, with 5Za adjusted 
in order to absorb the divergence. The UV finite self-energy then reads: 



nC(Q) = & tr / 7 At S (K) Y S (K - Q) + 5Z a \q 2 gT - Q»Q U \ ■ 



(42) 



The role of 5Za in the ring resummation is illustrated in Fig. El In general, the ring resummation 
does not account for all the next-to-leading order contributions. Indeed, the fermionic mass and 
field strength counterterms (5M and bZ^) are of order 1/Nf, and can be inserted in a fermion 
loop (see Fig. EJ) to generate a contribution of order N®. This contribution is needed to remove 
potential sub divergences as illustrated in Fig. [7| 



4 A more detailed analysis of the asymptotic behaviour can be found in the appendix of Ref. [3; 
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FIG. 5: Role of the countcrtcrm SZa- every time that a fermion bubble is inserted, one has to insert the 
corresponding counterterm 8Za to absorb the subdivergence. 



5M w 8z 




FIG. 6: Fermion "counterterms" eliminating next-to-leading order divergences in the bubble resummation 
of Fig.H 



Once all sub divergences are eliminated, there remains an overall, temperature independent 
divergence in the pressure, which can be eliminated by subtracting the vacuum pressure at next- 
to-leading order: 



P- P 



f 



NN f tT 



N n 



tr 



4 (SM + JfibZ^) S - [ (SM + fSZj,) S 
7k Jk e 

4 log^ + n) - / log (Gq- 1 + n vac ) 

Tq Jq 



(43) 



It is perhaps surprising that the first line of this formula did not appear in the general derivation 
given at the beginning of this section. This is because, the derivation there was made using bare 
quantities. If one rewrites Eq. (|38|) in terms of renormalized quantities, one generates these extra 
terms. Performing the Matsubara sums in Eq. (|43j) one obtains: 

P-Pl = NNf tr / {SM + fSZj,) {$ - M) a {K) - N g tr / e(q )n(\q \) Im log (Gq 1 + U) 

JK JQ 



N n 



tr 



log(Go 1 +n) - log (G^ + U^ c ) 



(44) 



The pressure in Eq. (|44l) contains a contribution carrying the counterterms, a second contribu- 
tion which is explicitly finite due to the presence of the factor n(|go|)> and a potentially divergent 




FIG. 7: Potential singularity with the structure of a fermion self-energy insertion. It appears when one of 
the fcrmionic momenta is kept fixed while all others momenta are taken to infinity. The structure of the 
diagram needed to absorb this potential singularity is exactly that of the diagram in Fig. 
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part which reads (A ~ B here means that the divergent parts of A and B coincide): 



(P-P c 



/\div 



tr 



log (g - 1 + n) - log (Go 1 + n vac ) 



(45) 



or, after performing the traces explicitly (see also |3J]), 

n tit h 



(P - Pi) div 



log 1 + 



-Q 2 + n v 



2 lo S I 1 4 



n 



e,th 



? 2 4n v 



(46) 



Using the asymptotic behaviour of Ht,th and H^th which can be found in the Appendix of Ref. 
one can show that the potential divergent term arises from the leading terms in the expansion of 
the logarithms: 



(P-P6 



/\div 



N„ 



th 4 2IT th 



2 J Qf -Q 2 + n v 



/ 1J th ' 



(47) 



or using the expression (|41[) for ITvi and Gf = N„/(2N): 



(P-P&) 



/\div 



NNf tr / (fi-M)ao(K) 
Jk 



g 2 C f / 7 I/ 5 (^-Q)7 M G f - c (Q) 



(48) 



As we shall see in Sect. II V Bl in the process of performing the analytic continuation in ko of the 
quantity in brackets in the formula above, one has to deform the contour of the Qe integration in 
order to avoid crossing singularities. However, here we are only interested in the UV contributions. 
Since contour deformation will only be necessary for Qe < 2/c = 2/co as we will see at the end of 
section UV B II we can identify the integral over Qe with the vacuum fermion self-energy (see Eq. 
(|69j) below), and write: 



(P - Pif v ~ NNf tr / {$- M) a (K) S b>vac [G vac ](^) . 

Jk 

We now add this divergent contribution to that of the counterterms 5M and 5Z^, and get 



(P-PJ) 



/\ct+di 



A 



NN f tr / (#-M)* (K) 5M + ^^ + S b , vac [G vac ](K) 



(49) 



(50) 



which can be made finite by a suitable adjustment of 5M and SZ^. 

In the massless case the divergences associated to the fermionic self-energies do not contribute 
to the pressure. As M vanishes, and = 0, it is clear that 5M does not contribute. One can 
verify that 5Z^ does not contribute in two ways. From (|50|) we have: 

(P-itf)* -^NNfSZ^tv [ f a Q {K)$ = 2nxNN f 8Z i) tr f f(\k \)5(K 2 )K 2 = 0. (51) 

Jk Jk 

Equivalently, one can rewrite the potentially divergent piece as: 



(P-P ( 



f\div 



NN 



f 



tr I fia (K)i: h>vac [G vac }(K) = 2irNN f [ f(\k \)5(K 2 )t vliC [G vllc ](K). 
Jk Jk 

(52) 

As we shall explain in section PV CI S vac [G vac ] = tr (^TS vac [G vac ]) vanishes on the light cone. 
This result is in one-to-one correspondence with the fact that, in the massless case, the large-Q 
contribution of the integrand in (|4*T|) , vanishes after integration over Q ji^, 13] • This is due to the 
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vanishing of the angular integral and requires that numerical computations use a cutoff respecting 
Euclidean symmetry. Although we shall eventually work with renormalized quantities obtained in 
the symmetry-preserving dimensional regularization, the numerical evaluations require a cutoff to 
eliminate the Landau pole of the large- Nf theory. As we shall discuss further in Sect. IIVI if this 
cutoff does not respect Euclidean rotation invariance, it leads to spurious contributions. 
With renormalized self energy II, Eq. (|44|) reduces in the massless case to @,Ei,E3: 



t N 
P-P/ = -ftr 



i log (G - 1 + n) - f log (Go 1 + n v 

To Jq„ 



(53) 



B. Large-iV/ limit of the self-consistent entropy 



The entropy could be obtained by taking a total derivative of the pressure with respect to the 
temperature. Using Eq. (|53|) for the pressure, this leads to 



d(P - Pp 
dT 



tr 



/ n(ft,)Im [G^ 



(54) 



where G 1 = G 1 + IT* and we subtracted the fermionic and the bosonic interaction-free contri- 
butions to the pressure Po — P$ + Pq- Note that the overall temperature- independent divergence 
that is present in the pressure disappears when considering the entropy. As we shall see, the only 
remaining divergences are associated with the self energies. This is, however, not manifest in Eq. 
(|54l) . where the second term is potentially divergent for ko — > — oo. 

Using that [G] is a two-loop diagram for which S' as defined in Eq. (|13j) vanishes identically, 
we have 



d(T$,[G\) 



dT 



where S is the 0(N7 1 ) fermionic self-energy obtained by calculating the one-loop diagram com- 
posed of a bare fermion propagator and a full large- Nf gauge boson propagator (see section ITVj) . 
On the other hand, 



tr 



G 



K 



^ 0) Ren*ImG + 2^ReSImS 



dT 



dT 



(55) 



d(T$*[G}) 



dT 



G 



d_ 
dT 

tr 



Tr(ILG) 



G 



d_ 

dT 



I n(£;o)Im(IUG) 

K 



G 



K 



dn(ko 
dT 



Im(n*G)+tr f n(ko)Im[G^ 
Jk V dT 



(56) 



Combining Eqs. and (jSS)) yields 

dn(k ) 



S-So 



tr 



-2tr 



K dT 
df(k , 



A 



dT 



Im log(G _1 G ) - ImnReG 
Re S Im , 



(57) 



where So denotes the ideal-gas limit of the entropy. 

This expression can in fact be derived also from the 2-loop-$-derivable expression for the en- 
tropy, Eq. (|14[). In the large- iVj limit, the full gauge boson self-energy II reduces to the fermion 
loop involving undressed propagators. Any fermion self energy diagram insertion (Fig. in II 
would bring in a factor g 2 without a factor Nf, leading to a subleading correction. However, one 
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fermion self-energy has to be included in tr log S 1-1 = tr log(S^" + E) in order to produce con- 
tributions of order Nj and N®. The integrand of the fermionic terms in the 2-loop entropy and 
density, Eqs. (|14j) and (|15|) . therefore simplifies according to 

Im log S- 1 - Im (S) Re (S) -► Im log S^ 1 + Im (ES* ) - Im E Re S + O^f 1 ) 

-> Im log Sq 1 + Re E Im 5 . (58) 

This enables us to recover Eq. (JBTj) for the entropy, while we get 

AA-M, = -2tr f ^4^ReSImS (59) 
Jk dn 

for the contribution of order N® to the density. 

In contrast to Eq. (|54|). Eq. (|57|) is a manifestly ultraviolet finite expression as soon as propaga- 
tors and self-energies are made finite through renormalization, and it naturally identifies bosonic 
and fermionic quasiparticle entropies, S = S b + Sf, with vanishing residual interactions S' . 



C. HTL approximation of the large- Nf entropy 

In the large- iVj limit, the two expressions for the entropy, Eqs. (|54|) and Eq. (|57|). are completely 
equivalent. However, this equivalence is broken once one starts approximating the gluon self- 
energy H* by its HTL limit which has different behaviour at large momentum. Doing so leads to 
uncancelled UV divergences in Eq. ([54)) . whereas the entropy formula ()57j) that was derived using 
the stationarity of the thermodynamic potential for self-consistent solutions remains finite. In the 
following we shall investigate the quantitative difference it makes to approximate the full large- Nf 
gluon self-energy by its HTL approximation and define 

Im log(G _1 G ) - ImftReG 

-2tr J ReE[G]W = S + S b + S fi (60) 

where II and G denote the gluon HTL self-energy and propagator, respectively. The numerical 
evaluation of this expression does not contain spurious divergences, because ReE[G] on the light 
cone is UV finite as we will show in section II V CI This is the advantage of using the self-consistent 
2PI formulation rather than a direct ring resummation. 

It is in fact instructive to evaluate the difference between Eq. (|6*0|) and a direct HTL approx- 
imation of d(P — Po)/dT as given by Eq. ()54j) . In the large- Nf limit the equivalence of the two 
entropy formulae stems from Eq. (|55|) and Eq. (|56|) which combine together into: 

tr ^^llmH,ReG-2tr J *!M ReEW = -tr ^n(fc )Im (g^\ . (61) 

In the HTL approximation, it is possible to derive a similar formula by exploiting the fact that 
Eq. 1)55)1 and Eq. (|56j) are valid for any G provided one maintains the functional relation of E to G. 
The value of n* does not need to be changed. Thus, if we replace G — ► G, E — > E[G], we obtain 

tr ^^^Imn*ReG-2tr ^ ^1 R e E[G] Im 5 = - tr ^n(^ )Im (g^) . (62) 



5n 



tr 



on{K ) 
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Now we add a common term to both sides in order to obtain the 2PI HTL entropy in the form: 



<5 = ^((P-Pi)lu^n 



dn(k ) 
dT 



Im(n-n*)ReG 



+ tr / n(ko) Im 
IK 



g ( <m fin* 



dT dT 



(63) 



The last two terms are the ones corresponding to the temperature dependent singularities in a 
direct HTL approximation on ring diagrams. 

Returning to the expression (J6(J|). we observe that the large- Nf limit leads to rather different 
simplifications in the bosonic and the fermionic contributions. In the bosonic part Sb, the NLO 
contributions <5II displayed in Fig. |2] do not appear because they are of order NJ X . The NLO 
contribution 5T, of Fig. |21 is also of order Ny 1 , but it does contribute to the fermionic entropy 
because there are Nf fermions. However, 5T, is not Dyson-resummed because for Nf — ► oo only 
one insertion of <5£ survives. 

Because of the factor Im Sq the fermion self-energy in (|57|) and Q6UJ) is evaluated on mass shell 
and because the integral is dominated by hard momenta, the fermionic contribution to the entropy 
can be identified with a weighted average over a momentum-dependent asymptotic quark mass. 
For the entropy 5 at zero chemical potential, we can define 



(S-S ( 



o)f 



-2tr 



df(k ) 



Re £ Im So 



-ANN 



dT 

df(k) Ml{k) 



NN f T 



dT 



2k 



6 



■M 



(64) 



although, as we shall see, M^ Q (k) and do not need to be positive. 

In perturbation theory we have in the massless case and at zero chemical potential 



N f Ml 



9cs 



+ 0(gi s ) 



CtT 1 



C, 



2N 



(65) 



where the calculation of the contribution oc g s requires HTL resummation as shown in Fig. |31 The 
latter is responsible for 3/4 of the plasmon term in the thermodynamic potential; the remaining 
1/4 of the plasmon term comes from the soft momentum regime of the bosonic entropy ()57|) . 

In the following we shall derive and evaluate numerically the quantities M^^k) and in the 
large- Nf limit. These will be compared with strictly perturbative results and HTL approximations 
where all the higher-order terms generated by HTL resummation are retained. This will allow us to 
test the proposal of complete HTL resummation not only with respect to the entropy of large- Nf 
QCD but also with respect to the asymptotic thermal quark masses. 



IV. FERMION SELF-ENERGY 



In this section we shall calculate the fermionic self-energy £ in the large- Nf limit as it is needed 
in the self-consistent entropy and number density formulae Eqs. Q57JI an d ()59|) . It corresponds to 
the one-loop diagram of Fig. El with an undressed fermion line and a dressed gluon propagator 
which resums the fermion bubbles in analogy to Fig. 0J The fermion self-energy requires mass and 



For the quark density a different weighted average would be relevant. 
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wave function renormalization: T,(K) = T,^(K) + SM + fiSZ^ with the bare self-energy given 
by 



V h {K) = -g 2 cd ^S (K-Q)YG^(Q). 
To 



It enters the fermionic contribution to the entropy in the second line of Eq. (|57|) as 



sf-sH 



-NN f 



1 

2Tk 



dT ReS(e fc ,A;) — — Re£(-e fc ,fc) 



with T>(K) from Eq. (|22j) . In the massless case this reduces to (see Appendix IB"]) : 



■NN 



I 



1 



Re S(fco = k,k) . 



(66) 



(67) 



(68) 



/k * 

We note here that the frequency carried by K in Eq. (|66jl is imaginary (Matsubara frequency). In 
contrast the evaluation of the entropy needs the discontinuity of the fermion self-energy across the 
real axis. The analytical continuation of (|66jl is described in detail in section II V Bl To that purpose, 
it is convenient to split the self-energy into a vacuum and a thermal piece £b = Sb, vac + ^th with 
respect to the overall loop momentum Q. The thermal piece is finite and its analytical continuation 
causes no troubles (see section TlV Aj) . In contrast, the correct evaluation of the vacuum part with 
respect to the integral over Q needs that we maintain Euclidean invariance (see section flVB|) . 
We thus compute it by replacing the discrete sum in (|6*6*)) by a continuous integral, leading to a 
4-dimensional Euclidean integral: 



X b ,vac W = ~9 2 Cf [ r S (K - Q) ^ G^(Q) . 



(69) 



This contribution also carries temperature dependences, but only through the gluon propagator. 
It contains UV divergences which in the massless case are washed out from the entropy formula 
as it is clear by looking at the contribution of the counterterms (once again the trace is performed 
over the Lorentz group): 



Re£ ct (±£ fc ,/c) = Retr 



{$ - M) (5M + $5Z^) 
AM5M . 



fco = ±£fe 



= ±4M 2 5Zj, - 4.M5M . (70) 

In section UV CI we give a direct proof of the finiteness of ReSb iV ac- In what follows, we drop the 
label "b" on the projected self-energy which reads: 



H(K) = tr [pE(K)] 



-9 2 C 



tr 



Q 



P = K-Q 



Evaluating the trace over 7 matrices one obtains 



£(*) 



-±9 2 C f 



Q 



K» P v + K v - K ■ P g»" A (P)G flu (Q) . 



(71) 



(72) 



In order to pursue the calculation, we need to specify a gauge. We shall use the Coulomb gauge 
from section Hi Bl By plugging G flu (Q) from Eq. (|TT)|) into Eq. (|72|). we obtain two components to 



t T {K) 



Ag 2 C 



-8g 2 C 



where &o ; Qo an d Po are evaluated at the Matsubara frequencies ko 
iuj — ioj n . 



fcopo + k-pj A (P) G L (Q), 

k p - (q • k) (p • q)l A (P) G T (Q) , (73) 

iu, qo = iuj n and po = iu m = 
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Matsubara sums and thermal contributions 



The contributions to X that we labelled as "thermal" are those which after performing the 
Matsubara sum in (|73|) contain a thermal factor n(|po|) or /(|po|) which vanishes when T — » 0. 
Note that these contributions do not account for the complete temperature dependences of S, since 
part of the latter is implicit in the dressed gluon propagator. The Matsubara sums are performed 
in Appendix IA 21 The result is: 



J L, th 



(K) 



9 2 Cf 



dppfp 



k 2 + p 2 , ( k— p 
1 H ~ log 



Akp 



k + p 



Ag 2 C f 



p 



kopo + k-p 



J T,th 



+ Vcv 

(if) = -8 5 2 C 
+ 8g 2 C f 



f 



ao (P) G^K-P) 
h (k - q ) + k • p A (A" - Q) <r L (Q) , 

fcoPo - (q • k)( P • q)] ffo(P) Gt(^ - P) 
fco (h ~ qo) - (q • k)(p ■ q)l A (K - Q) a T (Q) 



p 



(74) 



(75) 



witha (P) = e(p ) /| po | po(P) ando- L ,T(Q) = e(go) ™| ?0 | /°l,t(Q) and ^ = (iui,k). These integrals 
are UV finite thanks to the presence of the thermal factors. The continuation to the real axis 
ko = icoi — > ko = uj + ie and the light-cone limit uj — > k pose no problems, since the integrand 
is always well defined. After performing the angular integrals analytically, one can perform the 
remaining 2-dimensional integrals numerically. 

Since integrations in Eth are effectively cut off by the temperature, the necessity of introducing 
an ultraviolet cutoff below the scale of the Landau pole is no practical problem for temperatures 
T -C Al. The results will be independent on how such a cutoff is introduced. However, this will 
be different for the remaining "vacuum" contributions. 



B. Vacuum contributions and renormalization 

The remaining contributions S vac , which contain temperature dependences (only) through the 
spectral data of the gluon propagator, can be computed by replacing the discrete sum in lf75|) by 
a continuous integral, yielding a 4-dimensional Euclidean integral: 



k p + k • p 



A (P)G L (Q), 



£t c (*0 = ~8g 2 Cf / fcopo-(q-k)(p.q) A (P)G T (Q), 



(76) 



where K = {k ,k) = (iui,k), P = (po,p) = {iuJ m ,p), and Q = (qo,q) = {iuJ n ,q) have imaginary 
frequency to start with, and P + Q = K. Since our theory contains a Landau pole, we have 
to introduce a cutoff below that scale, and it is important that this cutoff is implemented in an 
Euclidean invariant way even when the result would be finite in dimensional regularization. This 
is not only important for quantities affected by the Landau pole but also for quantities free from 
the Landau singularity, as we illustrate in what follows. 
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1. Cutoff implementation example 



To give an example of difficulties encountered with a non Euclidean-invariant cutoff, consider 
the expression J{K) = I(K) - 1(0) with 

w = 4«F*r (77) 

This quantity appears in the calculation of T, vac [Go](K) and we would like to evaluate it on the 
light cone K 2 = 0. In dimensional regularization one can shift the integration momentum, and 
the result is [I(K) — 1(0)] — > [1(0) — 1(0)] = for arbitrary Euclidean momentum K. Hence, the 
light-cone value of J is zero. Is it possible to recover this result by using an explicit cutoff? 

Let us first introduce a cutoff that violates Euclidean invariance: a cutoff that only applies to 
three-dimensional momenta, not to the frequency (note that we assume imaginary values for qo 
and ko such that the integrand does not contain any poles along the integration path): 



(2vr) 4 7_ loo i J V-i '(qo-ko) 2 -(q 2 + k 2 -2qkcos9) 

= — ^A 2 + -^k 2 . (78) 
8vr 2 24vr 2 V ' 

The integrations over 8, qo, and q are straightforward (in this order), but the result we obtain, 
I(K) — 1(0) = fc 2 /24-7r 2 , is clearly in contradiction to the result of dimensional regularization. This 
result turns out to be independent of ko and can thus trivially be continued to Minkowski space 
leading to a non-zero value on the light cone. It turns out that this is a consequence of having 
introduced a cutoff procedure that violates Euclidean rotation invariance for the expression I(K). 

Thus, if one has to split frequency and momentum, one should still do it in a Euclidean invariant 
way: 

m = f A M f^ <* /' dx 



(2vr) 4 7 J-iy/J^ i J -i (qo ~ k ) 2 - (q 2 + k 2 - 

T lf A , dq {qo _ ko ) 2 -(q-k) 2 

d 11 , — - — - ~~T~ § 71 17Y2 TTTTvi- ( 79 ) 



(2ir)*kJ HH J_ iy ^Z^ i b (q - k ) 2 - (q + kf 

The angular integration along x = cos 9 poses no problem for imaginary qo = iu) n and ko = iui 
as the integrand is finite and continuous for all — 1 < x < 1. The integration along qo is more 
problematic as the integrand 

fa-fco) 2 -^ (8Q) 
(qo ~ k ) 2 -(q + k) 2 

has singularities at qo = ko ± \q ± k\ and branch cuts which might interfere with the integration 
path for an external Minkowski = (u>,k). The particular form of writing (|8U|) this way (i.e. not 
splitting up the logarithm \og(X/Y) into logX — logY) leads to an analytic structure in the qo 
complex plane as shown in Fig. |HJ a branch cut connecting the left two singularities and a branch 
cut connecting the rightmost two singularities. Note that if we start with purely imaginary ko, the 
(/o integration poses no problem since all singularities and branch cuts are away from the imaginary 
integration axis (except for q = k). 

In order to evaluate (|79|) on the light cone, we have to rotate ko from the imaginary to the 
real axis ko = ioJi — ► oj + ie. By doing so, all four singularities will move to the right, and one of 
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In* / A 2 ^2 
+Z V A — (J 


5o 


A B 


C D 




■AAA* 











FIG. 8: The branch cut structure for the qo integration in i|79[) with the logarithmic singularities at A = 
ko — \q + k\, B = ko — \q — k\, C = ko + \q — k\, D = ko + \q + k\. For purely imaginary ko (in the plot 
just between B and C), the integration path runs from —i^Jh 2 — q 2 to i\J h 2 — q 2 and does not cross any 
logarithmic branch cut. 




FIG. 9: As the rotation ko = iuj — > oj moves all the singularities at the same time, we should also deform 
the integration path in order to avoid crossing a singularity or branch cut. 

them will eventually cross the integration path when Re^o > \q — k\. We want I{K) to stay an 
analytic function. This is only possible if we avoid singularities and branch cuts by deforming the 
(numerical) integration path as shown in Fig. |UJ A deformation of a complex path that does not 
cross poles, singularities, or branch cuts, will not change the result of a complex integration. 

Note that q = k gives a "pinch singularity" , that is a point where in Fig |H] two singularities pinch 
the integration path. It turns out that in our case this singular point corresponds to a discontinuity 
of the first derivative of the q integrand, but is otherwise harmless for the q integration. Also 
note that when q > 2k = 2ko, we do not have to deform the integration path as the point B = 
ko~\q~k\ < in Fig. I^Jstays on the left side of the integration path even on the light cone. We can 
choose the integration path such that all path deformation is performed for Qe = \/~Qo + Q 2 < ^i 
and no path deformation needed for Qe > Ai as long as the intermediate cutoff Ai > 2k. 
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-iy/A 2 - q 2 



qo 



A 



-z^A 2 - g 2 



FIG. 10: A general propagator with self-energy II will introduce additional cuts and poles along the real qo 
axis. In this case, the integration has to go through the origin and ko has to be shifted slightly off the real 
axis to ko = u> + ie. 

Performing the numerical integration with the modified path indeed gives the desired result 
I(K) — 1(0) = also on the light cone, in accordance with dimensional regularization. Blindly 
integrating across the branch cut without deforming the integration path would give a wrong, 
non-vanishing result on the light cone. 

2. Full cutoff implementation 
We can now perform the angular integration in the vacuum contributions (|76|) and obtain 



-4g 2 C f 



2vr 



dq 



V^ 3 ? dq 



Gl(Q) 



{qo - k )(qo - 3&o) -q 2 + k 2 



4kq 



log 



\K) 



-8g 2 C f 



2tt 



A 



dq q 



o 

2\2 



(qo-ko) 2 -(q-k) 2 
(qo-k ) 2 -(q + k) 2 

■ KoY + 



^ Gt(Q) 



{(qo ~ k ) 2 -k 2 Y- q 2 (Ako(k - q ) + q 2 ) , 
+ 8kq* bg 



(go - kp) 2 + q 2 
2q 2 

(qo ~ kp) 2 -(q- k) 2 
(qo-ko) 2 -(q + k) 2 



k 2 



(81) 



The bare propagator Ao(P) = Ao(K — Q) gives the same logarithmic expression of the form 1)80(1 as 
in our example. But Gl(Q) and Gt(Q), which contain the gluon self-energies III and IIt, provide 
additional branch cut structures along the real qo axis such that the path deformation is further 
restricted as shown in Fig. ^3 The integration path goes from qo : A 2 — q 2 — > — ► ko(= 
LO + ia) — > lyj A 2 — q 2 . 

Another numerical subtlety is involved: Keeping e small but fixed will limit the upper integration 
bound of the qo integration to i\J A 2 — q 2 > ie (see Fig. 110(1 which means that the q- integration is 
limited to < q < VA 2 - e 2 . The error introduced by this turns out to be at least of the order 
Q(e 2 ). 
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C. Renormalization 

1. Large-Nf limit 

As already mentioned, the disappearance of the counterterms 8Z^ and SM from the entropy 
formula indicates that Re£ vac has to be finite on the light cone. Since the propagator G used to 
compute this quantity contains a Landau pole, we have to be more specific about what we call 
"finite". Re 

S V ac is computed with an explicit Euclidean invariant cutoff A <§C Al. The finiteness 
of Re S vac means that it is insensitive to A in a broad region of momenta T< A<C Al. In order 
to understand this more explicitly, we can use the linearity of Re S V ac with respect to G in order 
to write Re S vac = Re S vac [G vac ] + Re S vac [G — G vac ]. Using the asymptotic behaviour of II t h, one 
can check that the second term of this expression is insensitive to A. Thus we only need to check 
the finiteness of ReE vac [G vac ] on the light cone. 

As explained in the previous section, the use of an Euclidean invariant cutoff allows one to obtain 
the same results that one would obtain in dimensional regularization if there were no Landau pole. 
If we suppose that there is no Landau pole, then it is possible to show that ReE vac [G vac ] = 
on the light cone. Indeed this is a gauge invariant quantity and one can compute it in any 
gauge. Choosing a covariant gauge, the cancellation is straightforward according to covariance and 
dimensional regularization. Of course we cannot rigorously use dimensional regularization because 
of the Landau pole. However using an explicit Euclidean invariant cutoff, it is possible to check 
that ReE vac [G vac ] = on the light cone. 

To do so, we start from equations (|81|) and transform the integration variables to 4-dimensional 
variables according to 



dq q 



dq 



f(qo,q) 



A f'TV 

dQE Qe / d£ sin 2 £ f(iQE cos£, Q E sin£) . 
o Jo 



(82) 



The integrand can then be expanded for large Qe and, since in the free case G^{Qe) = G^{Qe) = 
—Gl(Qe)/ sin 2 £ = (Q 2 E + ^^(Q^))" 1 is independent of £, we can perform the £ integration and 
obtain 
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(83) 



The first term of each expansion is logarithmically divergent if transverse and longitudinal parts 
are taken separately. Only when adding the components together to form Y, va " c (K) 
E™ c (.ff), the logarithmic divergence vanishes on the light cone k^ = k 2 : 



Re£ vac (K) 



-Ag 2 C f 



2tt 2 



dQE QeG YSiC {Qe) 
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Ql 
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15Q| 



01 i 



(84) 
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2. HTL approximation 

Unlike the case of large Nf, the disappearance of the counterterms 5Z^ and 5M from the 
entropy formula cannot be taken as an indication of the finiteness of ReS vac [G]. This is because 
G is built out of approximations of Feynman diagrams whose asymptotic behaviour is changed 
with respect to the large- Nf case. If one writes ReS vac [G] = ReS vac [Go] + ReT, v3iC [G — Go] one 
is able to check that the first term of this expression is zero on the light cone (using an Euclidean 
invariant cutoff). For the second term, since the asymptotic behaviour of II is different from that 
of the large- Nf result n t h, we need a direct check. In principle one can perform this check as in 
the previous section, just Eq. (|83|) would have to be changed when using HTL propagators, since 
Gt(Qe, £,) and Gi(Qe, £) now also depend on the angle £. One can easily convince oneself that the 
angular dependence does not affect the leading logarithmic divergence though, as both transverse 
and longitudinal HTL propagators can be expanded in the form 

^ = wA l+m t +0 (if)}- (85) 

The additional angular dependence / (£) is suppressed by a factor m 2 D /Q 2 E at large Qe and therefore 
does not influence the leading logarithmic divergent behaviour of Eqs. Q83|) and ()84j) : the leading 
divergence still vanishes on the light cone. 



V. NUMERICAL RESULTS 
A. Numerical implementations 

The evaluation of equation (j7fifl requires three consecutive numerical integrations (of which one 
is given by the complex path described earlier) in order to obtain one value of £(u; = k; g 2 s )- With 
our implementation this takes of the order of 10 hours on a current PC (3GHz). This has to be 
integrated over to obtain the fermionic contribution to the entropy SWmion(fl , cff) ^ or one P°i n t °f 
g 2 s . Part of the code was therefore ported to run in parallel and final results could be obtained on 
the ECT* Teraflop Cluster within a couple of weeks that would have taken years on a standard 
PC. 

It turns out that the numeric cancellation in the integrand (|76jl for large q and (imaginary) 
qo works best if one integrates over 4-dimensional Euclidean spheres at fixed Qe = V~Qo + I 2 
first. No path deformation is needed as long as Qe > 2k as was discussed at the end of section 
IIV B II Unfortunately, the path deformation for 4-dimensional angular variables turns out to be 
much more involved (the corresponding deformed path could become infinitely long due to inverse 
trigonometric functions), such that a hybrid approach seems a good compromise between simple 
path-deformation for small Qe and good convergence properties for larger Qe- 

A Mi r A 

4„ _ / j4„ , / .4. 



JO J\i 



Al 7 2 H^~ q2 dq r 



J-i^/XJ-c 



2 I 



dn 2 + / dQ E Q% / dn 3 



J dn 3 (86) 



where Ai is an intermediate cutoff to be chosen such that 2k < Ai <C A (in our implementation 
we chose Ai = 2.1k), and J d&2 an d J dQs denote the 3- and 4-dimensional angular integrations. 
Practically, yet another cutoff Ai < A2 < A is introduced, above which the angular integration 
j dil.3 is performed analytically on the high-temperature series expansion. This series expansion 
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FIG. 11: Asymptotic thermal quark mass squared (real part of the fermionic self-energy on the light-cone) 
as a function of k/T for g^ s (irT) = 4. The exact large- Nf result is compared to the NLO-HTL calculation 
and its relevant average value. The renormalization scale /ij^g is varied around the FAC-m scale by factors 
of 2 (see text). 

is also used to show that neither large- Nf nor HTL gluon self-energies contain logarithmically 
divergent pieces for high momentum shells. 

In order to correctly integrate all peaks along the complex path, the routine for complex path 
integration is enhanced by providing information about analytically or numerically known positions 
in the complex plane of poles of the propagator, o>l(<z) and cjt(<z) |4o |. singularities of the logarithms 
and other singular points. The integration path is divided into smaller segments in the vicinity of 
such points where the integrand can rapidly change its value by orders of magnitude. 



Figures 111! 1121 and EH show the results of a numerical calculation of the asymptotic thermal 
quark mass squared M^(k) = 2kReT,(ko = k) for three different values of the coupling g 2 s (/^ms = 
ttT) = 4, 9, and 16, and normalized by T 2 Cf/Nf. (Recall that this quantity is of order N7 , and 
it contributes to the entropy only because there are Nf fermions.) The exact (nonperturbative) 
result obtained in the large- Nf limit is given by the full lines. 

As anticipated, the asymptotic thermal quark mass squared is a nontrivial function of momen- 
tum. It is maximal for k ~ T and decays for large k 3> T. The fact that it vanishes and even 
becomes negative at very small momentum is not to be taken seriously. M^k) is referred to as 
asymptotic thermal quark mass, because only at hard k >T does this have the interpretation of a 
quasiparticle mass; for smaller k one would have to search for a self-consistently determined pole 
of the dressed quark propagator rather than evaluate the quark self-energy on the tree-level mass 
shell, i.e. the light-cone. M^ Q (k) however also becomes negative at very large momenta, k S> T. 
This means that the dispersion law for hard fermionic modes eventually turns space-like, which 
does not necessarily signal an inconsistency of the theory. In fact, the group velocity remains 
smaller than the speed of light for all k > T. Negative M^ a (k) however opens up the possibility 
for Cerenkov radiation and thus means that at the corresponding values of k there is a qualitative 
difference to the leading-order result = ^g 2 s CfT 2 /Nf. 



B. Asymptotic thermal quark masses 
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Ml(k)N f /(T*C f ) 




FIG. 12: Same as Fig.[TT]for g% s (irT) = 9. 
Ml(k)N f /(T 2 C f ) 




FIG. 13: Same as Fig.QUfor 9es(^ T ) = 16 . 

The next-to- leading order result as obtained in HTL-resummed perturbation theory (i.e., re- 
placing the full gauge boson propagator in the quark self-energy by its HTL approximation) is 
given by the dashed lines in Figs. I11H13I Qualitatively, the HTL result is similar to the exact 
large-iVy result. Quantitatively, like every perturbative result, it has a dependence on the renor- 
malization scale. The HTL result is displayed for a range of renormalization scales /%fg centered at 
/2FAC-m = exp(i — 7e)7tT ps 0.9267rT which is the scale of fastest apparent convergence with respect 
to the effective mass parameter rn-E in dimensional reduction jj| ■ The upper (dotted) limiting lines 
of the shaded bands correspond to \ x /zfac-iih while the lower (dashed) limiting lines correspond 
to 2 x /xpAC-m- The renormalization scale dependence is seen to increase with the coupling g^ s , 
and also in the domain where k S> T. For large values of k the HTL result for the asymptotic 
mass squared turns more quickly negative than the full result, but in fact at large values of k the 
modes are severely suppressed by the Bose-Einstein distribution factor. The range k ~ T is more 
relevant physically (in particular for the calculation of the entropy to be discussed presently), and 
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FIG. 14: Comparison of the averaged asymptotic thermal quark mass squared, M^, in various approxima- 
tions as explained in the text. The renormalization scale p-j^g is varied around the FAC-m scale by factors 
of 2. 



here the agreement with the exact large- iVj result is much better. 

For the computation of the fermionic contribution to the entropy, the asymptotic thermal mass 
appears in form of the weighted average given in Eq. (|64jh This averaged asymptotic thermal 
mass is given for the HTL approximation by the flat horizontal bands labelled (NLO-HTL). 
For the entire range of couplings that we can consider without becoming sensitive to the scale 
of the Landau pole, the averaged asymptotic quark masses are compared in Fig. 1141 The full 
line is again the exact large- iVy result, given as a function of g 2 ff (7rT). The short-dashed lines 
correspond to the leading-order HTL value with renormalization scale varied as before. The 
perturbative next-to-leading order result (|65|) is given by the longer-dashed lines. Finally the full 
NLO-HTL-resummed result, not truncated at order g^ s , is given by the dash-dotted lines. 

As one can see in Fig. 1141 the NLO-HTL-resummed result represents a considerable improvement 
over the HTL-resummed result truncated at order g^ s for g 2 ff > 4, which corresponds torao/T > 1. 



C. Numerical results for the entropy and discussion 

Figure ^] shows the numerical results of the entropy calculation. The full line is the entropy 
density S = (dP/dT)^ as it has been obtained earlier 0, |34], EH by a numerical derivative of 
the pressure P from Eq. (|43|) for massless fermions m = 0. It is a very non-trivial numerical test 
that the result for the entropy as obtained from the ^-derivable two-loop approximation, Eq. (|14|). 
reproduces exactly the same result in the large- Nj limit, as has been discussed in section IHI Al 
Indeed, only after implementing the correct path deformation for the calculation of S, we have 
been able to reproduce the entropy to three or four digits, where the accuracy was limited only 
by the high calculation cost for the multi-dimensional integrals. In Figure El the two curves lie 
perfectly on top of each other for the whole range of couplings displayed. 

We show the results in a ran ge o f couplings g 2 ff (7rT) < 36 where the influence of the Landau 
pole can be neglected. As in 

0,13,0 we obtain differences in the entropy on the percent level 
for the largest couplings ffc ff (/%jg = vrT) « 36 by varying the numerical cutoff A 2 = aA 2 in the 
range a = 1/4 .. 1/2. 
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FIG. 15: Entropy in the large-A^ limit separated into a bosonic and a fermionic part, corresponding re- 
spectively to the first and second line of Eq. (|57|l . The NLO-HTL approximation to the two parts of the 
entropy is depicted in dashed lines, where the renormalization scale pr^g is varied around the FAC-m scale 
by factors of 2. The combined NLO-HTL result shows remarkable agreement with the exact large- Nf result 
for all couplings with HTL Dcbye mass rho < 2.5T. 
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FIG. 16: Entropy in the largc-iV/ limit, comparing the exact large-TV/ result to the strictly perturbative 
expansion through order g^ s , the NLA result, and the NLO-HTL result. The renormalization scale is 
varied around the FAC-m scale by factors of 2. 

According to Eq. (|57|) , the large- Nf entropy as well as the HTL approximation thereof is com- 
posed of a bosonic and a fermionic contribution. The fermionic one is entirely given by the weighted 
average (|64|) of the asymptotic thermal quark mass that we have discussed above and it is repro- 
duced in Fig. El by the line marked "fermion" , together with its NLO-HTL approximation. The 
calculation of the bosonic contribution is less demanding computationally. The result is shown by 
the line marked "boson" in Fig. 1151 The HTL approximation to the bosonic contribution does not 
require NLO corrections to the asymptotic bosonic masses in the large-A^y limit, and it has been 
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calculated completely already in Refs. [21]. |2J|, Evaluated for the range of renormalization scales 
considered here it gives the band below the line marked "boson" . Perhaps fortuitously, the errors of 
the bosonic contribution to the HTL approximation of the entropy are opposite in sign from those 
of the fermionic contribution. The sum total turns out to reproduce the exact result with astound- 
ing accuracy up to g^ s ~ 16. For larger coupling the renormalization scale dependence quickly 
becomes enormous, however the central value determined by the optimized scale fij^ = /iFAC-m 
remains amazingly close to the exact result for all values of g c R. In Ref. j3f| a similarly successful 
approximation was constructed by using optimized renormalization scales on a perturbative result, 
which however required to include all contributions up to and including order g® s . 

In Fig. the exact and the NLO-HTL result is also compared to simpler approximations. 
The short-dashed lines give the strictly perturbative result up to and including order <^ ff , the 
longer-dashed lines labelled "NLA" (next-to-leading approximation) show the result of comparing 
the bosonic HTL entropy with a perturbative approximation to the averaged asymptotic thermal 
quark mass presented before in Ref. [3(|. The latter corresponds to the large- iVj limit of the 



HTL-resummed entropies obtained in Refs. |2ll . |22J, |2j]. While this does lead to an important 



improvement for couplings such that rh^/T < 1.5, it is remarkable that a complete HTL resum- 
mation in the asymptotic thermal mass can push the range where HTL resummation works well 
up to ni£)/T ~ 2.5. 



VI. CONCLUSIONS AND OUTLOOK 

We have calculated the large-Aj limit of the entropy of ultrarelativistic gauge theories by 
evaluating separately the contributions from bosonic and fermionic quasiparticles. Since in the 
large-iVj limit, the 2-loop ^-derivable approximation becomes exact, this allowed us to assess the 
error that is being made when fully dressed propagators and self-energies are replaced by their 
HTL approximations. 

In the case of the bosonic contributions, these approximations correspond to replacing a full 
fermion loop in the dressed gluon propagator by its HTL approximation. The fermionic contri- 
butions to the entropy, on the other hand, involve the momentum-dependent asymptotic thermal 
quark mass where the HTL approximation has to be carried to next-to-leading order. In both 
cases the HTL approximation turned out to give remarkably good results when compared with 
the exact large- Nf results, even for fairly large coupling. Combining these contributions, the final 
result for the entropy in the NLO-HTL approximation turned out to be amazingly accurate up 
to rho/T ~ 2.5, and even beyond when the renormalization scale is fixed by the requirement of 
fastest apparent convergence of the electric mass in dimensional reduction. Its quality is then 
comparable to optimized perturbative results including terms through order g® s . This seems to 
support the conclusions of Ref. 4^, H3] where optimised dimensional-reduction results for QCD to 



order g 6 log(g) were found to agree well with available lattice data and in turn with the estimates 
from the HTL entropy. 

The improvement achieved by including the full momentum dependence of asymptotic thermal 
masses and keeping all effects of the resummation of hard thermal loops is encouraging for further 
developments of this approach. A straightforward extension of the present calculations would be the 
inclusion of fini te q uark chemical potential. The calculation of quark densities and quark number 
susceptibilities |5ll . I5H ] for small- Nf QCD within NLO-HTL resummation in fact requires 
exclusively the asymptotic quark thermal masses that were calculated in the present work for zero 
chemical potential. By also calculating the full momentum dependent asymptotic gluon masses, 
one could finally complete the existing HTL results 2l], 2^, [i^] for the entropy of full QCD in the 
2-loop ^-derivable approximation. 
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APPENDIX A: MATSUBARA SUMS 
1. Gluon self-energy 

The gluon self-energy to leading order in Nf is given by: 

K(Q) =9eS tr^ ^S (K)YS (P), (Al) 

with P = K — Q. It is convenient to split this self-energy into a vacuum and a thermal piece lib = 
lib, vac + n t h (the label "b" stands for "bare" and is here to remind that the vacuum piece contains 
UV divergences). The vacuum piece is simply obtained by replacing the discrete Matsubara sum 
by a continuous integral, thus leading to a 4-dimensional Euclidean integral: 

n r v ac(Q)=5e 2 ff tr f rS (K)YS (P). (A2) 
JKts 

The thermal piece is obtained after performing the Matsubara sums and separating the thermal 
dependent pieces. A convenient way to proceed is to introduce the spectral representation (J7J) for 
each of the fermionic propagators: 

nf(Q) = s c 2 ff tr III Y($-M)p (K)Y(r-M)p (P)^n — -^7 

(A3) 

where uj m = lo\ — io n and qo = iu n . One can now easily perform the Matsubara sum and extract 
the thermal dependent part: 

ly 1 -/(fc ) + /(p ) _ e(fco)/(|fc |)-e(po)/(|po|) | ^ 

P i (h - iuJi){p - iui + iu n ) po-k + iLu n p -k + iuj n 

where we have used f(ko) = 8(—ko) — e(ko) f (\ko\) . The thermal part of the gluon self-energy now 
reads (we do not use the label "b" since this contribution is UV finite): 

nST«»=Af / / / /y-M)Mq/(f-M)Mq ' (W(IMI :: W/fal) . (At) 

Jpo Jk Jk Po-k + qo 

Using the spectral representation for the fermion propagator, the parity properties of IP" together 
with a straightforward change of variables (which is justified since the cut-off can be sent to infinity 
in the thermal contribution), one obtains: 

D£(Q) = 2 5c 2 ff tr f <f{$-M)a Q {K)i'S {K-Q), (A6) 
Jk 

with a (K)=e(k )f(\k \) Po (K). 
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2. Fermion self-energy 



In this section, we compute the fermionic Matsubara sums: 



E L (K) = -Ag 2 C f 



Q 



k po + k • p 



MP) Gl(q) 



S t (K) = Sg'CfT fc Po-(q-k)(p-q) A (P)G T (Q) 



(A7) 



with P = K — Q. To that aim, we use the spectral representations (0) and (|18|) . The longitudinal 
part gives: 



t L (K) = Ag 2 C f 



4g 2 C f 



q J po 
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(A8) 



where iu m = itoi — iu> n and ko = iuj\. The Matsubara sums give: 



Iv I 

P „ po - m + iw„ 



n-po+iw, = /-po = e(po)/(bo|) + vac , 
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n qo + /— p 



e(goMkol) ~ e(ffo)/(bol) 
Po + 9o - i^i 



+ vac , 



(A9) 



where we have used f(p ) = 0(-po) ~ e(Po)/(|Po|) and n(q ) = -9(-q ) + e(q )n(\q \) . Then, we 
have: 



S L ,th(*Q = 4< 7 2 C / 



Ag 2 C } 
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q 2 



Po (P)e(po)f(\ Po \) 



q J p a J q 
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(A10) 



We can now use the spectral representations backwards to write (we change variables in the second 
line q — > p, this is possible since the cut-off can be sent to infinity in the thermal contributions): 



kp 



p q 



MP) 



ko po + k • p 



a (P) G L (K-P) 



+ 4g 2 C f / k (ko - q ) + k • p A (K - Q) a L (Q) , 



(All) 



with (To(-P) = e(po) /|p | Po(P) and o"l(Q) = e(<7o) n |g | Pl(Q)- There is a static contribution which 
can be computed almost analytically since it does not depend on the gluon propagator: 

f k • p Q 2 Cf f°° f 1 kpx 

WCf Jp -^ MP) = —Jo dPPU jj X W+P^2kp-x 



g 2 C~ 



S 



7T- 



dppfp 



k 2 + p 2 . ( k — p 
1 H n log 



Akp 



k + p 



(A12) 
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In the same way, one obtains for the transverse piece: 




fc Po-(q-k)(p-q) a (P)G T (K-P) 



+ 8g 2 C f 



I 

JQ. 



ko (ko - q ) - (q • k)(p • q) A (K - Q) a T (Q) 



(A13) 



with cr T (Q) = e(g ) n M pr(Q). 



APPENDIX B: RELATION BETWEEN £+ AND £_ 



and S_ from Eq. (|24|) can be related to each other by application of the Schwarz reflection 
principle: It can be shown that b(ko = iu>i,k) = b*(iui,k) is a purely real quantity in Euclidean 
space (i.e. along the imaginary axis ko = iui), as is ia(iu)i, k) = —ia*(iuji, k). The Schwarz reflection 
principle implies that in a rotation from Euclidean to Minkowski space ko = iu>i — > uj + ie on the 
one hand and ko = iu>i — > —co + ie on the other, the values of b can be related via b(ko = uj + ie,k) = 
b*(—iv + ie, k). The same holds for ia{ui + ie, k) = —ia*(—io + ie, k). Combining the two results as 
in S-t = b ± a (and E± = b* ± a*) gives 



This result is of course also valid for finite e, i.e. all uj > and e > 0. Particularly for e — > this 
implies 
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